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ABSTRACT 

We use a high resolution iV-body simulation to study the velocity bias of dark matter halos, 
the difference in the velocity fields of dark matter and halos, in a flat low-density ACDM model. 
The high force, 2/i~ 1 kpc, and mass, 10 9 /i _1 Mq, resolution allows dark matter halos to survive in 
very dense environments of groups and clusters making it possible to use halos as galaxy tracers. 
We find that the velocity bias &„,12 measured as a ratio of pairwise velocities of the halos to 
that of the dark matter evolves with time and depends on scale. At high redshifts (z ~ 5) halos 
move generally faster than the dark matter almost on all scales: b v ^n(r) « 1.2, r > 0.5/i _1 Mpc. 
At later moments the bias decreases and gets below unity on scales less than r w 5/i Mpc: 
b v ,i2(r) ~ (0.6 — 0.8) at z = 0. We find that the evolution of the pairwise velocity bias follows 
and probably is defined by the spatial antibias of the dark matter halos at small scales. One-point 
velocity bias b v , defined as the ratio of the rms velocities of halos and dark matter, provides a 
more direct measure of the difference in velocities because it is less sensitive to the spatial bias. 
We analyze b v in clusters of galaxies and find that halos are "hotter" than the dark matter: 
b v = (1.2 — 1.3) for r — (0.2 — 0.8)r„j r , where r vir is the virial radius. At larger radii, b v decreases 
and approaches unity at r = (1 — 2)r vir . We argue that dynamical friction may be responsible 
for this small positive velocity bias (b v > 1) found in the central parts of clusters. We do not 
find significant systematic difference in the velocity anisotropy of halos and the dark matter. The 
dark matter the velocity anisotropy can be approximated as (3(x) — 0.15 + 2x/(x 2 + 4), where 
distance x is measured in units of the virial radius. 

Subject headings: cosmology:theory - large-scale structure of universe - methods: numerical 



1. Introduction 

Peculiar velocities of galaxies arise due to the 
gravitational pull of surrounding overdense regions 
and therefore reflect the underlying density field. 
The statistical study of galaxy velocities is im- 
portant in cosmology since it can be used as a 
tool to constrain cosmological models. The con- 
nection between theoretical predictions and the 
observed statistics usually requires an additional 
quantity: the difference between galaxy and dark 
matter velocities, termed the velocity bias. The 
situation with predictions of the velocity bias is 



rather confusing. There is a wide range of esti- 
mates of the velocity bias. Values change from 
strong antibias with galaxies moving twice slower 
than the dark matter (Gelb & Bertschinger 1994; 
Klypin et al. 1993), to almost no bias (Klypin et 
al. 1998 (KGKK); Ghigna et al. 1998), to slight 
positive bias (Diaferio et al. 1998; Okamoto & 
Habe 1999). Following Carlberg (1994) and Sum- 
mers, Davis, & Evrard (1995) we distinguish two 
forms of the velocity bias. The one-point velocity 
bias b v is defined as the ratio of the rms velocity 
of galaxies or galactic tracers to that of the dark 



1 



matter: 



b v = 



Cgal 

ctdm ' 



(1) 



where the rms velocity a is estimated on some 
scale. Traditionally, this measure of the velocity 
bias is used for clusters of galaxies. Two-particle 
or pairwise velocity bias 6^12 compares the rela- 
tive velocity dispersion in pairs of objects sepa- 
rated by distance r: 



b v ,\2 — 



E (r) 



Cdm,dmM 



(2) 



The pairwise velocity dispersion (PVD) was of- 
ten used to complement the analysis of the two- 
point spatial correlation function. At small scales, 
the cosmic virial theorem Peebles (1980) pre- 
dicts that the PVD of galaxies should be pro- 
portional to the product of the mean density of 
the universe and the two-point correlation func- 
tion. The PVD of galaxies has been estimated for 
the CfA catalog (Davis & Peebles 1983; Zurek et 
al.1994; Somerville, Davis & Primack 1997) and 
recently for the Las Campanas Redshift Survey 
by Landy, Szalay, & Broadhurst (1998) and Jing, 
Mo & Borner (1998). The latter two studies gave 
363±44km/s and 570±80km/s, respectively, for a 
l/i _1 Mpc separation. Jing & Borner (1998) show 
that the discrepancy between these two studies is 
due to the difference in treatment of the infall ve- 
locities. The value of er gig as well as the infall ve- 
locities depend on which regions (clusters or field) 
are included in the surveyed sample. 

The PVD of the dark matter, <7dm,dm, has also 
been estimated for a variety of cosmological mod- 
els (e.g., Davis et al. 1985; Carlberg & Couch- 
man 1989; Carlberg, Couchman & Thomas 1990; 
Klypin et al.1993; Colin, Carlberg, & Couchman 
1997; Jenkins et al. 1998). If galaxies were 
a random sample of the mass distribution, we 
would expect that <r g . g were approximately equal 
to <7dm,dm- Davis et al. (1985) showed that in this 
case an ilo = 1 model with erg = 1 produces a 
PVD that is too large compared to observations. 
Here os is the rms of mass fluctuation estimated 
with the top-hat window of radius 8/i _1 Mpc. This 
is an example of a model which needs some kind 
of bias to be compatible with the observations. 

The notion of the pairwise velocity bias b v _i 2 
was introduced by Carlberg & Couchman (1989). 
They found that the dark matter had a PVD a 



factor of two higher than that of the simulated 
"galaxies" . In a further analysis, Carlberg, Couch- 
man, & Thomas (1990) suggested that an 17o = 1 
model with as — 1 could be made consistent with 
the available data for b v ^2 ~ 0.5 (velocity an- 
tibias) and almost no spatial bias. Estimates of 
the pairwise velocity bias are in the range of 0.5- 
0.8 (Couchman & Carlberg 1992; Cen & Ostriker 
1992; Gelb & Bertschinger 1994; Evrard, Sum- 
mers & Davis 1994; Colin, Carlberg, & Couchman 
1997; Kauffmann et al. 1998a, b). Differences be- 
tween the estimates (especially the early ones) can 
be attributed to some extent to numerical effects 
( "ovcrmcrging problem" ) and to different methods 
of identifying galaxy tracers. Only recently N- 
body simulations achieved a high dynamic range 
in a relatively large region of the universe neces- 
sary for a large number of galaxy-size halos to sur- 
vive in clusters and groups (e.g.,KGKK, Ghigna et 
al. 1998; Colin et al. 1999). The estimates of the 
pairwise velocity bias start showing a tendency for 
convergence. For example, results of Kauffmann 
et al. (1998a, 1998b) for a low-density model with 
a cosmological constant and results presented in 
this paper for the same cosmological model agree 
reasonably well in spite of the fact that we use very 
different methods. Results point systematically to 
a antibias b Vi \2 = 0.6 — 0.7. 

One-point velocity bias for clusters and groups 
of galaxies tells a different story. Values of b v 
are typically larger than those for b Vi \2 and range 
from 0.7 to 1.1 (Carlberg & Dubinski 1991; Katz 
& White 1993; Carlberg 1994; Ghigna et al. 
1998; Frenk et al. 1996; Metzler & Evrard 1997; 
Okamoto & Habe 1999; Diaferio et al. 1998). 
Carlberg & Dubinski (1991) suggested that if the 
pairwise velocity antibias is significant, galaxies in 
clusters should have orbital velocities lower than 
the dark matter. However, this may not neces- 
sarily be true. In this paper (see also, for ex- 
ample, Kauffmann et al. 1998a) we argue that 
galaxy tracers do not need to move slower in clus- 
ters to have the pairwise velocity bias b Vt i2 < 1. 
In particular, we find that while 6^12 < 1 for halos 
in our study, the halos in many clusters actually 
move somewhat faster than dark matter. Ghigna 
et al. (1998) also do not detect a significant differ- 
ence between the orbits of DM particles and halos. 
They find that the cluster radial velocity disper- 
sion of halos is within a few percent of that of 
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the DM particles. Okamoto & Habe (1999) used 
hundreds of galaxy-size halos in their simulated 
cluster. They are able to compute the halo veloc- 
ity dispersion profile. Their results suggest that 
in the range 0.3 Mpc ^ r ^ 0.6 Mpc halos have 
a velocity dispersion slightly larger than that of 
the DM particles. Diaferio et al. (1998) using a 
technique that combines iV-body simulations and 
semi-analytic hierarchical galaxy formation mod- 
eling also find that galaxies in clusters have higher 
orbital velocities than the underlying dark matter 
field. They suggest that this effect is due to the 
infall velocities of blue galaxies. We find in this 
paper a similar effect: galaxy-size halos are "hot- 
ter" than the dark matter in clusters. 

The paper is organized as follows. In § 2 brief 
descriptions of the model, simulation, and group- 
finding algorithm are given. In § 3 the DM and 
halo PVDs as well as the corresponding velocity 
bias are computed at four epochs. We take a sam- 
ple of the most massive clusters in our simulation 
and compute an average halo and DM velocity dis- 
persion profile. A cluster velocity bias is then de- 
fined and computed. A discussion of the main 
results are presented in § 4. The conclusions are 
given in § 5. 

2. Model, simulation, halo-finding algo- 
rithm 

We use a flat low-density model (ACDM) with 
flo = 1 — 51a = 0.3 and ag = 1. Cluster mass 
estimates (e.g., Carlberg et al. 1996), evolu- 
tion of the abundance of galaxy clusters (e.g., Eke 
et al. 1998), baryon fraction in clusters (e.g., 
Evrard 1997), and the galaxy tracer two-point 
correlation function(e.g., Colin et al. 1999; Ben- 
son et al. 1999) favor a low-density universe with 
fio ~ 0.3 (see also Roos & Harun-or-Rashid 1998). 
On the other hand, various observational deter- 
minations of h (the Hubble constant in units of 
100 km sec -1 Mpc" 1 ) converge to values between 
0.6-0.7. Our model was set to h = 0.7 which gives 
an age for the universe of 13.4 Gyr in close agree- 
ment with the oldest globular cluster age deter- 
minations (Chaboyer 1998). The approximation 
for the power spectrum is that given by (Klypin & 
Holtzman 1997). The adopted normalization of 
the power spectrum is consistent with both the 
COBE observations and observed abundance of 



galaxy clusters. 

The Adaptive Refinement Tree code (ART; 
Kravtsov, Klypin & Khokhlov 1997) was used to 
run the simulation, as described by Colin et al. 
(1999). The simulation followed the evolution of 
256 3 dark matter particles in a 60/i _1 Mpc box 
which gives particle mass of 1.1 x lO 9 ft. _1 M0. The 
peak force resolution reached in the simulation is 
~ 2/i~ 1 kpc. The mass resolution is sufficient for 
resolving and identifying galaxy-size halos with at 
least 30 particles. The force resolution allows ha- 
los to survive within regions of very high density 
(as those found in groups and clusters of galax- 
ies) . In dense environment of clusters the mass of 
halos is not well defined. Therefore, we use the 
maximum circular velocity 

y^J^A) 1 ", (3) 

V / max 

where M(< r) is the mass of the halo inside radius 
r, as a "proxy" for mass. 

Halos begin to form at very early epochs. For 
example, at z ~ 6 we identify > 3, 000 halos with 
maximum circular velocity, V max , greater than 
90 km/s. The numbers of halos that we find at 
z = 3, 1, and are 14102, 14513, 10020, respec- 
tively. We use a limit of 90 km/s on the circular 
velocity which is slightly lower than the complete- 
ness limit - (110 - 120)km/s (Colin ct al. 1999) 
of our halo catalog. This T^nax value increases the 
number of halos quite substantially (a factor of 
two as compared with the limit of 120 km/s), and, 
thus, reduces the statistical noise. We checked 
that our main results are only slightly affected by 
partial incompleteness of the sample. 

Our halo identification algorithm, the Bound 
Density Maxima (BDM; see KGKK), is described 
in detail elsewhere (Klypin & Holtzman 1997). 
The main idea of the BDM algorithm is to find 
positions of local maxima in the density field 
smoothed at the scale of interest (20/i _1 kpc). 
BDM applies physically motivated criteria to test 
whether a group of DM particles is a gravitation- 
ally bound halo. The major virtue of the algo- 
rithm is that it is capable of finding both isolated 
halos and halos orbiting within larger dense sys- 
tems. Cluster-size halos were also located by the 
BDM algorithm. The physical properties of a sam- 
ple of the twelve most massive groups and clus- 
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Table 1 

Physical Parameters of Clusters' 1 



M v { r CTZD Vmax Rvir f^halo 



h^Me 


(km/sec) 


(km/sec) 




Knax > 90km/; 


6.5 x 10 14 


1645 


1402 


1.43 


246 


2.4 x 10 14 


1022 


910 


1.28 


132 


1.9 x 10 14 


992 


831 


1.17 


98 


1.6 x 10 14 


975 


789 


1.11 


95 


1.4 x 10 14 


887 


747 


1.05 


58 


1.3 x 10 14 


887 


730 


1.02 


55 


1.1 x 10 14 


831 


695 


0.98 


45 


1.0 x 10 14 


820 


680 


0.95 


33 


9.9 x 10 13 


789 


673 


0.94 


74 


9.7 x 10 13 


789 


668 


0.94 


60 


9.3 x 10 13 


753 


659 


0.92 


67 


8.3 x 10 13 


720 


635 


0.89 


64 



a Column description: (1) virial mass of the cluster; (2) 3D velocity 
dispersion of dark matter particles; (3) maximum circular velocity; (4) 
cluster radius; (5) number of galaxy-size halos with V max > 90km/s. 



ters 1 are shown in Table 1. The total number of 
clusters chosen for the sample is a compromise be- 
tween taking a relatively large number of clusters, 
so that we could talk about cluster average prop- 
erties, and using clusters with a relatively high 
number of halos. This cluster sample is used to 
compute the average DM and halo velocity disper- 
sion profiles as well as the average DM and halo 
velocity anisotropy profiles. 

3. Results 

3.1. The pairwise velocity bias 

The three-dimensional pairwise velocity disper- 
sion PVD is defined as 

a 3 2 D (r) = (v 12 2 )-(v 12 ) 2 (4) 

where V12 is the relative velocity vector of a pair 
of objects separated by a distance r and brack- 



The cluster number 8, in descending order in mass, was 
excluded from the sample because it has a group close to 
it that produces too much disturbance to the cluster. 



ets indicate averaging over all pairs with the sep- 
aration r. Figure 1 shows the PVD for the dark 
matter, a^D,dm, at four epochs (top panel). At 
1 /i _1 Mpc the radial PVD is about 1100 km/s at 
z = 0. For the same cosmological model Jenkins 
ct al. (1998) find a radial PVD of - 910 km/s. 
Jenkins et al. used slightly lower normalization 
for the model (cr 8 = 0.9) and used a bigger sim- 
ulation box (L box = 141.3/i _1 Mpc). When the 
differences in as are taken into account the Jenk- 
ins et al. value increases to 1,120 km/s. Thus, 
both estimates roughly agree. 

The ratio of the halo and the dark matter 
PVDs, the pairwise velocity bias 6„,i2, is shown 
in the bottom panel of Figure 1. All halos with 
Vmax > 90 km/s were included in the computa- 
tion. At very early epochs and on large scales ha- 
los tend to move faster than the dark matter. At 
later moments the pairwise velocity bias becomes 
smaller than unity. It is interesting to compare 
the evolution of b Vy i2 with the changes in the spa- 
tial bias for the same cosmological model (Colin 
ct al. 1999). The spatial bias is defined as the 
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Fig. 1. — Top panel: Three-dimensional pairwise rms 
velocity of the dark matter at four different epochs in- 
dicated in the Figure. Bottom panel: Pairwise veloc- 
ity bias for galaxy-size halos with circular velocities 
V m ax > 90km/s. Curves are labeled in the same way 
as in the top panel. At very early epochs and on large 
scales halos tend to move faster than the dark matter. 
At later moments the pairwise velocity bias becomes 
smaller than unity. 

square root of the ratio of correlation functions 
[Chh(^) / 'CdmM] 1 ^ 2 • I n general, the biases evolve in 
the same way. At high redshifts both biases are 
positive (6 > 1) and decline as the redshift de- 
creases. At low redshifts biases dive below unity 
(antibias) and stop evolving. In spite of similari- 
ties, there are some differences. The pairwise ve- 
locity bias becomes less than unity at z = 3 on 
scales below 3/i _1 Mpc. At the same redshift the 
spatial bias is still positive on all scales. 

Colin et al. (1999) and Kravtsov & Klypin 
(1999) interpret the evolution of the spatial bias 
as the result of several competing effects. Statisti- 
cal bias (higher peaks are more clustered) tends to 
produce large positive bias and explain bias evo- 
lution at high redshifts. At later epochs halos of 
a given mass or circular velocity become less rare 
and start merging inside forming groups of galax- 
ies. Both effects lead to a decrease of bias. The 
merging becomes less important as clusters with 
large velocity dispersions form at z < 1. This 
results in a very slow evolution of the halo corre- 
lation function and bias. It is likely that the same 



processes define the evolution of the pairwise ve- 
locity bias. The differences can be explained by 
the known fact that the PVD is strongly dom- 
inated by few largest objects (e.g., Zurek et al. 
1994; Somerville, Davis & Primack 1997): merg- 
ing of halos inside forming groups at z — 3 results 
in fewer pairs with large relative velocities and in 
velocity antibias on w l/i _1 Mpc scales. If this in- 
terpretation is correct, the pairwise velocity bias 
mostly measures the spatial bias, not the differ- 
ences in velocities. 

3.2. The velocity anisotropy (3 

A sample of 12 groups and clusters (see Table 
1) was used to compute various average cluster ve- 
locity statistics. In order to reduce the noise in the 
profiles because of the small number of clusters in 
the sample, we double the sample by using also the 
same clusters at slightly different time z = 0.01. 
For each cluster the halo distances to the clus- 
ter center are divided by the corresponding cluster 
virial radius (normalized distances). The halo ve- 
locities (averaged in spherical bins) are divided by 
the corresponding cluster circular velocity at the 
virial radius (normalized velocities). In Figure 2 
we show radial profiles, in normalized units, for 
halos and DM: the mean radial velocity (v r ), the 
radial (oy) and the tangential (cr t ) velocity dis- 
persions. All halos are given equal weight. We 
have accounted for the Hubble flow when we com- 
pute <j r and er t (so, proper, not peculiar velocities 
are used); no correction for the mean radial ve- 
locity was made. The trend in both the velocity 
dispersion and the anisotropy velocity is slightly 
affected if the mean radial velocity is subtracted 
at distances << 0.6 and it is not affected at all at 
smaller distances. 

The velocity anisotropy function 

(3 = 1- <j\ l2oi (5) 

is presented in the bottom panel of Figure 2 for 
halos and for DM. For pure radial orbits (3=1, 
while an isotropic velocity dispersion implies (3 = 
0. The two lines added to the panel show a fitting 
formula (Carlberg et al. 1997). 

4r 

(3 = /3 m ^— -+0o (6) 

for two pairs of parameters (/3 TO , (3q): (0.65,0.) 
and (0.5,0.15). The first set of parameters gives 
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Fig. 2. — Velocity profiles averaged over all clusters in 
Table 1 excluding cluster #8. Curves are for the dark 
matter and different symbols are for halos with circu- 
lar velocity Vmax > 90km/s. Distances to the clus- 
ter centers are divided by the corresponding cluster 
virial radius and halo velocities (averaged in spherical 
bins) are divided by the corresponding cluster circu- 
lar velocity at the virial radius. Errorbars show 1- 
sigma errors of the mean. Top panel: Mean radial 
velocity (dot-dashed curve), the radial velocity disper- 
sion (full curve), and the tangential velocity dispersion 
(dashed curve) of the dark matter. Open squares, 
filled, and open circles show the radial velocity, the 
radial velocity dispersion, and the tangential velocity 
dispersion respectively for halos. Bottom panel: Ve- 
locity anisotropy for halos (open triangles) and for the 
DM (solid line) . The dot-dashed and three-dot-dashed 
lines represent the fitting [3 — 4r(3 m /(r 2 + 4) + (3o for 



two pairs of (/3„ 
tively. 



/3o): (0.65,0) and (0.5,0.15), respec- 



a better approximation for halos. It explicitly as- 
sumes that j3 = at the center. The second set of 
parameters allows a small anisotropy at the cen- 
ter. It provides a better fit for the dark matter. 
Note that while the halos have a tendency for more 
isotropic velocities (with possible exception of the 
center), the difference between halos and the dark 



matter is not statistically significant. 

The variances of ay and at are computed using 
standard expressions for errors; for example, for 

M4 



var(a r 



f4 



(7) 



where fi 2 = J2i( v r,i ~ v r ,i) 2 and ^ 4 = J2i( v r,i ~ 
Mr,*) j and n is the number of halos. The statistical 
error is, thus, given by the square root of var(a r ). 
The variance of /3 is given by 



[var{(3)f 



var(a 2 t ) 
2ol 



var(a 2 ) 
2ai 



(8) 



3.3. The cluster velocity bias 



The three-dimensional velocity dispersions for 
both halos and DM are shown in the top panel of 
Figure 3. The bottom panel shows the cluster ve- 
locity bias, defined here as b v = o- 3D j la i /a 3D , dm . 
It is surprising that halos in clusters appear to 
have larger, by about 20%, velocity dispersions 
than the DM particles (positive bias). The trend 
is the same regardless of what type of velocity dis- 
persion (3D, tangential or radial) we use in the 
velocity bias definition. There is almost no bias in 
the very center of clusters. However, the b v value 
of the innermost bin increases if we exclude the 
"cD" halos (defined as those halos which lie within 
the inner ~ 100ft - : kpc radius and have maximum 
circular velocities greater than about 300 km/s). 
Their exclusion increases the positive velocity bias 
to 1.22, a value which is comparable to that found 
in the adjacent bin. 

The cluster positive velocity bias is robust to 
changes in the limit of the circular velocity V max . 
Only the innermost bin experiences significant 
changes when this limit is increased. For example, 
when we increase Vmax from 90 km/s to 150 km/s 
(more massive halos are chosen) the value of b v in 
the innermost bin reduces to 0.6. This favors a pic- 
ture in which the central regions of clusters large 
galaxy-size halos feel more the slowing effect of 
the dynamical friction. All the other bins (within 
the virial radius) continue to show small positive 
velocity biases. The positive velocity bias is also 
robust to changes in the number of clusters of the 
sample. For instance, one might suspect, that the 
most massive cluster weights so much that it al- 
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Fig. 3. — Top panel: 3D velocity dispersion profiles 
for halos (filled circles) and DM (solid line) in units of 
the mean virial velocity. Bottom panel: Cluster veloc- 
ity bias profile. Errors correspond to 1-sigma errors of 
the mean. 

ters the statistics 2 . This is true to some degree, 
but is does not change the "sign" of the bias. For 
example, when we exclude this cluster and take 
Knax = 150 km/s, all bins continue to show pos- 
itive bias (within the virial radius) except the in- 
nermost bin where b v = 0.5. The results for the in- 
nermost bin should be taken with caution because 
the effects of overmerging may still be present in 
the central 100/i _1 kpc of the clusters. 

The difference in velocity dispersions of halos 
and dark matter particles indicates that their ve- 
locity distribution functions (VDF) arc different. 
We have examined both differential and cumu- 
lative VDFs for the analyzed clusters and found 
that the halo VDFs are generally skewed towards 
higher velocities as compared to the dark matter 
VDF, at r/r mr < 0.8. The two VDFs are ap- 
proximately the same for larger radii. The ob- 
served differences in the velocity distribution may 
be caused either by the differences in velocity fields 
of infalling halos and dark matter (if, for exam- 
ple, halos are accreted preferentially along fila- 
ments resulting in orbits of higher ellipticity) or 
by effects of dynamical friction operating on ha- 
los, but not on dark matter, in clusters. The dy- 



2 In fact, the most massive cluster of our simulation has had 
a recent major merger and halos may still have large ("over- 
heated") velocities (e.g., Katz & White 1993) 



namical friction may affect the slowest halos more 
efficiently because the dynamical friction time is 
proportional to the cube of the halo velocity. The 
slowest halos may therefore merge more efficiently 
thereby skewing the velocity distribution of the 
surviving halos towards higher velocities. 

One could ask whether or not this positive clus- 
ter velocity bias persists in the next set of twelve 
clusters or groups, in descending order in mass 
(with virial masses below those clusters shown in 
Table 1). Because this new sample of clusters have 
an average virial mass lower than the average mass 
of clusters of Table 1, dynamical friction is ex- 
pected to operate more efficiently (e.g., West & 
Richstone 1988; Diaferio et al. 1998). The num- 
ber of halos per cluster or group in this new sample 
is small, we therefore use the whole group velocity 
dispersion. We find integral b v values which are 
in general lower than one, and in some cases there 
are groups that exhibit a strong velocity antibias 
(ratios close to 0.6). This is contrary to what we 
find for the clusters of Table 1, where the majority 
of clusters have an integral positive velocity bias. 

4. Discussion 

Literature on the velocity bias is very extensive 
and results are often contradictory. In this sec- 
tion we review some of the published results and 
compare them with our results. There are some 
reasons for the chaotic state of the field. One of 
them is the confusion of two different notions of 
the velocity bias - the single-point b v and the pair- 
wise b v ,i2 biases. The biases have different nature, 
and, thus, give different results. Another source 
of confusion is the way how galaxies are identi- 
fied or approximated in theoretical models. When 
we combine this uncertainty with many physical 
processes, which we believe can create and change 
velocity bias, the situation becomes rather com- 
plicated. 

Velocity profiles seems to be the easiest part 
of the picture. In this paper we present results, 
which are less noisy and are based on a more ho- 
mogeneous set of clusters than in most of previous 
publications. Our results on the average cluster 
profiles for the dark matter (v r and ay) roughly 
agree with the results of Cole & Lacey (1996), 
Tormcn, Bouchet, & White (1997), and Thomas 
et al. (1998). For example, Tormen, Bouchet, 
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& White (1997) find a DM velocity anisotropy 
Pdm ^ 0.2 at r/r V i r 0.1 and (3dm — 0.5 at 
r/r V i r ~ 1, which is close to our results. The 
structure of galaxy clusters in various cosmologies 
is analyzed in detail by Thomas et al. (1998). 
From a total sample of 208 clusters they choose 
a subsample which shows no significant substruc- 
ture. They find a more isotropic averaged beta 
profile (l3 d m ~ 0.3 at r/r 180 = 1) in their ACDM 
model. The differences between our result and 
theirs can be accounted for the fact that their clus- 
ters were selected not to have significant substruc- 
ture. More substructure in a cluster likely means a 
more anisotropic cluster. The (3 value at the clus- 
ter center (innermost bins) is around 0.1, which is 
close to our results. 

Pairwise velocity bias is very sensitive to the 
number of pairs found in rich clusters of galax- 
ies. Removing few pairs may substantially change 
the bias. Thus, it mostly measures the spatial 
bias (or antibias) and is less sensitive to real dif- 
ferences in velocities. The value of 6^12 that we 
find at z — is typically higher than previous es- 
timates reported in the literature, computed for 
the Oo = 1 CDM model (Carlberg & Couchman 
1989; Carlberg, Couchman & Thomas 1990; Gelb 
& Bertschinger 1994; Summers, Davis, & Evrard 
1995) Some of the results are difficult to compare 
because the pairwise velocity bias is expected to 
evolve with time and vary from model to model. 

The first interesting result of this paper, that 
comes out from the evaluation of 6^12 at very high 
redshift, is that the halo PVD can be greater than 
that of the DM. This positive velocity bias had 
not been detected before (but see below) partly 
because of the lack of simulations with very high 
resolution that could overcome the overmerging 
problem. This result is surprising in part because 
halos are expected to be born dynamically cool 
3 . In fact, this is one of the reasons given in the 
literature to explain the present-day pairwise ve- 
locity bias (e.g., Evrard, Summers & Davis 1994). 
The other is the dynamical friction (e.g., Carlberg, 
Couchman & Thomas 1990). We offer the fol- 
lowing explanation to this positive velocity bias. 
Those halos that are formed at very high redshift 
come from very high density peaks. They are dy- 

3 Halos tend to form near the peaks of the DM density dis- 
tribution, (e.g., Frenk et al. 1988). 



namically cooler than an average DM particle from 
the region where they were born in, but hotter 
than most of the matter. The pairwise velocity 
bias b Vt i2 rapidly becomes smaller than one at non- 
linear scales. As time goes on, the mergers inside 
forming groups reduce the number of high veloc- 
ity halos, while velocities of DM particles increase. 
As the result the average halo random relative ve- 
locities are reduced below that of the DM. 

Using a semi-analytical method to track the 
formation of galaxies Kauffmann et al. (1998a; 
1998b) also find a pairwise velocity bias greater 
than one at high redshifts. They find that the 
galaxy PVD is greater than that of the DM at 
z > 1.1 (their figure 11, tCDM model). A 
b v , 12 > 1 is expected at higher redshift in their 
ACDM model as well. 

Single-point velocity bias appears to be the 
most difficult and controversial quantity. It is im- 
portant because it is a more direct measure of the 
velocity differences. It still depends on the spatial 
bias, but to much lesser degree as compared with 
the pairwise bias. An interesting result was found 
when we evaluated the average cluster halo veloc- 
ity dispersion profile and compared it with that 
of the DM particles: within the virial radius halos 
move faster than the dark matter. 

We believe that the explanation for this fact 
comes from a combination of two known physi- 
cal mechanisms: the dynamical friction and the 
merging of halos. One may naively expect that 
the dynamical friction should always slow down 
halos, which must result in halos moving slower 
than the dark matter particles. This is not true. 
While on a short time-scale the dynamical friction 
reduces velocity of a halo, the halo may decrease 
or increase its velocity depending on the distribu- 
tion of mass in the cluster and on the trajectory 
of the halo. For example, if a halo moves on a 
circular orbit inside a cluster with the Navarro- 
Frcnk- White profile, its velocity will first increase 
as it spirals from the virial radius to 2.2R S , where 
R s « (200 — 300) kpc is the characteristic radius of 
the core of the cluster. The halo velocity will then 
decrease at smaller radii. When the halo comes 
close to the center of the cluster it merges with 
the central cD halo, which will have a tendency 
to increase the average velocity of remaining ha- 
los. It appears that the Jeans equation provides a 
better tool for understanding the velocity bias. 
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We will use the Jeans equation as a guide 
through the jungle of contradictory results. It can- 
not be used more than a hint because it assumes 
that a cluster is stationary and spherical, which is 
generally not the case. If a system is in a station- 
ary state and is spherically symmetric, the mass 
M(< r) inside radius r is related to the radial ve- 
locity dispersion ay: 

M(<r)= r 4A, (9) 

* = -(Tng + Irg + W), (io) 

where p is the (number) density profile, and (i is 
the velocity anisotropy function. The left-hand- 
side of this equation (the total mass) is the same 
for both halos and the dark matter. Thus, if the 
term A is the same for the dark matter and halos, 
then there should be no velocity bias: halos and 
the dark matter must have the same ay. Numer- 
ical estimates of the term A are inevitably noisy 
because we have to differentiate noisy data. Nev- 
ertheless, we find that the value of the term A for 
halos is systematically smaller than for the dark 
matter. This gives a tendency for ay to be larger 
for halos. In turn, this produces a positive veloc- 
ity bias. The main contribution comes typically 
from the the logarithmic slope of the density: the 
halo density profile is shallower in the central part 
as compared with that of the dark matter. The 
halo profile is shallower likely because of merg- 
ing in the central part of the cluster, which gave 
rise to a central cD halo found in each of our clus- 
ters. We note that while the Jeans equation shows 
the correct tendency for the bias, it fails to repro- 
duce correct magnitude of the effect: variations of 
the term A are smaller than the measured velocity 
bias. 

One can also use the Jeans equation in a dif- 
ferent way - as an estimator of mass. We have 
computed M(< r) for our average cluster using 
both DM and halos. At (r/r vir ) = 0.25, where b v 
is close to its maximum, the halo mass determina- 
tion is larger than that of the DM by a factor of 
1.4. This is due to the larger halo velocity disper- 
sion. Because the term A is actually higher for DM 
by about 10%, the overestimation is reduced from 
1.56 to 1.4. As the distance to the cluster center 
approaches the virial radius the mass overestima- 
tion disappears. At the virial radius both mass 
estimations agree, essentially because (5, oy, and 



the sum of the logarithm derivatives are the same 
for both halos and DM, and are within (10 — 15)% 
of the true mass. 

Using the Jeans equation for a spherically sym- 
metric system and assuming an isotropic velocity 
field, Carlberg (1994) showed that a cool tracer 
population, b v < 1, moving inside a cluster with 
a power-law density profile (the density profile for 
the tracer is also assumed to be a power-law), pro- 
duced a mass segregation. That is, the tracer pop- 
ulation had a steeper density profile. We can in- 
vert this reasoning and say that a more centrally 
concentrated halo distribution produces a velocity 
antibias. We do not find this kind of mass segre- 
gation in our halo cluster distribution. In fact, we 
see the opposite - halos are less concentrated than 
DM. Dynamical friction along with merging pro- 
duces a lack of halos in the center of the cluster. 
This very likely explains differences between our 
and Carlbcrg's results for the velocity bias. 

Carlberg & Dubinski (1991) simulated a spher- 
ical region of 10 Mpc radius and 64 3 DM parti- 
cles. They were unable to find galaxy-size ha- 
los inside cluster at z — because of insuffi- 
cient resolution: softening length was 15 kpc in- 
stead of ~ 2/i~ 1 kpc needed for survival of halos 
(KGKK). Their identification of "galaxies" with 
those DM particles which were inside high-density 
groups found at high redshift, may have produced 
a spurious cluster velocity antibias. Using differ- 
ent galaxy tracers Carlberg (1994) also found an 
integral cluster velocity bias lower than one. This 
result could still be affected by numerical resolu- 
tion (e = 9.7/i _1 kpc). Evrard, Summers, & Davis 
(1994) run a two- fluid simulation in a small box, 
Lbox — 16 Mpc, and stopped it at z = 1. Each DM 
particle had a mass 9.7 x 10 8 M Q and an effective 
resolution of 13 kpc (at z — 1). The initial condi- 
tions were constrained to assure that a poor clus- 
ter could form in their simulation. Their "globs" 
(galaxy like objects) exhibit a lower than one ve- 
locity bias. This velocity bias appeared not to 
depend on epoch and mass. Their velocity antib- 
ias qualitatively agrees with our results for groups 
and poor clusters. At the same time, their value 
for the pairwise velocity bias agrees with our re- 
sults. 

Metzler & Evrard (1997) use an ensemble of 
two-fluid simulations to compute the structure of 
clusters. Unfortunately, their simulations do not 
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have high mass resolution to allow the gas in their 
simulations to cool and form "galaxies" (and then 
they could also allow for some feedback). In- 
stead, they use a high-density peak recipe to con- 
vert groups of gas particles into galaxies particles. 
They find a one-point "galaxy" velocity bias that 
depends on cluster mass: the higher the cluster 
mass is the higher the b v value. We find a simi- 
lar result when we do the analysis of the velocity 
bias cluster by cluster 4 . Their ensemble-averaged 
bias parameter is 0.84. Their recipe for galaxy for- 
mation produces a galaxy number density profile 
which is steeper than that of the DM. This is likely 
the reason why they find a b v value lower than one 
(Carlberg 1994, see above). 

Frenk et al. (1996) simulated a Coma-like clus- 
ter with a P 3 M + SPH code that includes the 
effects of radiative cooling. The mass per gas 
particle is 2.4 x 10 9 M Q with a softening param- 
eter e — 35 kpc of the Plummer potential. Their 
galaxies have two extreme representations: one 
as a pure gas clumps and the other as lumps of 
the stellar component. They find a mass segrega- 
tion in both representations - galaxies are more 
clustered than DM toward the center of the clus- 
ter which is not seen in our halo distribution 5 . 
Once again, according to Carlberg (1994) analy- 
sis, this would result in a one-point velocity bias 
lower than one (b v ~ 0.7). Because of a strong 
cooling, their "galaxies" can acquire high density 
contrasts, which helps galaxies to survive inside 
cluster. At the same time, poor force resolution 
(35 kpc) could have affected their results. 

There are two studies where b v values greater 
than one are obtained. Okamoto & Habe (1999) 
simulate a spherical region of 30 Mpc radius us- 
ing a constrained random field method. They use 
a multi-mass initial condition to reach high res- 
olution. Their high-resolution simulated region, 
where the cluster ends up, has a softening length 
e = 5 kpc and mass per particle m ~ 10 9 M Q . 
They find a cluster velocity bias lower than one 
only in the innermost part of the cluster where dy- 
namical friction is expected to be more efficient. A 
small positive bias (b v > 1) is found in the range 
0.3 Mpc < r < 0.6 Mpc. Based on the previ- 

4 On individual clusters we take only integral velocity dis- 
persions 

5 The reader might want to compare the Fig. 11 in Frenk et 
al. (1996) with the Fig. 2 in Coli'n et al. (1999) 



ous work by Kauffman et al. (1998a), Diafcrio 
et al. (1998) study properties of galaxy groups 
and clusters. They also find that galaxies in clus- 
ters are "hotter" than the underlying dark matter 
field. They suggest that this effect is due to the 
infall velocities of blue galaxies. Infall could ex- 
plain the positive velocity bias of the outermost 
bin (within the virial radius) of our Figure 3, but 
it definitely cannot account for the b v > 1 value 
seen in the inner bins (the mean radial velocity is 
close to zero for both DM and halos in the three 
innermost bins). 

There are several differences between our sim- 
ulation and those mentioned above. First, some 
of the papers cited above simulate only a region 
which ends up as a cluster. So, they have struc- 
ture for only one cluster. The single-cluster one- 
point velocity bias could not represent an average 
velocity bias, found using a sufficiently large sam- 
ple of clusters. For example, if our small positive 
velocity bias is influenced by non-equilibrium clus- 
ter features, then when one selects a cluster which 
is in good dynamical equilibrium (this could be 
defined, for example, by the absence of substruc- 
ture in the cluster) and computes the one-point 
velocity bias, it could be biased toward low values 
(b v < 1) because dynamical friction have had more 
time to operate. Second, we simulate a relatively 
large random volume that gives us many clusters 
in which effects such as tidal torques, infall, and 
mergers are included naturally. A cluster simu- 
lated region or a random large region but without 
enough resolution may not have a sufficiently large 
number of galaxy tracers and, thus, introduce high 
statistical errors. Our relatively large number of 
halos in clusters reduces significantly the statisti- 
cal errors in the computation of b v and makes them 
suitable to the determination of, for example, the 
radial dependence of the velocity bias. Third, in 
view of the Okamoto & Habe (1999) and Ghigna 
et al. (1998) results, and our own results, it seems 
that numerical resolution not only plays an impor- 
tant role in determining the whole cluster velocity 
bias value (both spatial and velocity bias inter- 
vene to affect its value) but it is also important in 
determining the radial dependence of b v (almost 
pure velocity bias). 

What could account for the small positive ve- 
locity bias that we see in our average cluster? We 
have examined both the differential and the cumu- 
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lative radial velocity distribution functions. We 
use the radial velocity to highlight any contribu- 
tion of infall velocities to the velocity bias. The 
cumulative radial velocity distribution function is 
shown in the Figure 4 for four different radial bins. 
In the top-right panel (the innermost bin) we see 
a higher fraction of low-velocity halos at small v r 
values. This is due to central cD halos, which 
move very slowly relative to clusters themselves. 
At large v r values we observe the contrary - a 
higher slope, which means that there are many 
fast moving halos. If we do not include the cD 
halos, the velocity bias becomes larger then unity 
even in the central radial bin. However, as we no- 
ticed earlier in section 3.3, a velocity antibias can 
appear in the central bin, if the value of Vmax is 
increased. It is clear that the deficiency of low and 
moderate v r halos produces the positive velocity 
bias measured at r = (0.2 — 0.8)r„j r (see the top- 
left and the bottom-right panels). We have used 
the Kolmogorov-Smirnov test in order to evaluate 
whether or not the halo and the DM velocity dis- 
tribution functions are statistically different. We 
find that the probability that these functions were 
drawn from the same distribution is smaller than 
0.01 in all radial bins that are within the virial 
radius. As mentioned in § 3.3, the dynamical fric- 
tion may have affected the slow moving halos more 
significantly because the dynamical friction time- 
scale is proportional to the cube of the halo ve- 
locity. It is thus expected that low-velocity halos 
merge sooner than their high-speed counterparts, 
thereby skewing the VDF toward high- velocity ha- 
los 6 . Infall could also be an important source of 
positive velocity bias for the outermost bins. 

5. Summary 

1. We have found that galaxy-size halos have a 
time- and scale- dependent pairwise velocity bias. 
At high redshifts (z ~ 5) this bias is larger than 
unity (w 1.2). It declines with time and becomes 
f» 0.6 — 0.8 at z — 0. The evolution of the pairwise 
velocity bias follows and probably is defined by 
the spatial bias of the dark matter halos. These 
results are in qualitative agreement with those by 

e It should be also kept in mind that as halos move to or- 
bits of smaller radii they could acquire higher velocities 
because the DM velocity dispersion increases toward the 
cluster center 
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Fig. 4. — Cumulative radial velocity distribution 
functions in four radial bins for halos (dotted line) and 
DM (solid line). 

Kauffmann et al. (1998b). 

2. We have evaluated the velocity anisotropy 
function (3(r) for both halos and DM particles. 
For both halos and DM (3 is a function that in- 
creases with radius and reaches a value of ~ 0.5 at 
the virial radius. The difference between this value 
and that found by Thomas et al. (1998) likely can 
be explained by the fact that Thomas et al. (1998) 
selected a sample of clusters which had little sub- 
structure. Our simulations indicate that the halo 
velocity anisotropy closely follows (but lies slightly 
below) that of the underlying dark matter. 

3. Halos in our clusters move faster than DM 
particles: b v = (1.2 — 1.3) for r — (0.2 — 0.8)r vir . 
This result disagrees with many previous esti- 
mates of the cluster velocity bias. This difference 
appears to be due to differences in numerical res- 
olution. More work needs to be done to settle the 
issue. Nevertheless, it is encouraging that Diafcrio 
et al. (1998) and Okamoto & Habe (1999) found 
results similar to ours. 

4. Usual argument that dynamical friction 
slows down galaxies and, thus, must produce ve- 
locity antibias is not correct. Galaxy tracers in 
clusters move through an environment which has 
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a steep density gradient. A sinking halo may cither 
increase or decrease its velocity depending on the 
distribution of cluster mass and on the trajectory 
of the halo. A combination of the dynamical fric- 
tion and merging appears as the most compelling 
hypothesis which could account for our small pos- 
itive velocity bias. 
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